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I.  INTRODUCTION 


Small  laboratory  devices  such  as  vented  chambers^or  blowout  guns 
have  a  long  history  in  understanding  gun  barrel  wear.  Compared  with 
large-caliber  gun  tests,  such  devices  need  less  propellant,  need  less 
room  for  safe  firing,  and  have  inserts  which  can  be  removed  for  post¬ 
firing  analysis.  An  unanswered  question  is  how  to  relate  wear  in  the 
devices  with  wear  in  large-caliber  guns.  A  recent  .JANNAF  workshop  con¬ 
cluded  that  laboratory  devices  cannot  evaluate  wear-reducing  additives, 2 
and  that  test  devices  in  three  laboratories  have  given  contradictory 
results  with  the  same  propellants. 3,4 

This  report  attempts  to  relate  erosivity  from  propellants  in  the 
BRL  37mm  blowout  gun  to  wear  expected  in  large  caliber  guns.  Empirical 
formulas5,6  presume  the  wear  is  exponentially  proportional  to  the  maxi¬ 
mum  bore  surface  temperature.  Bore  surface  temperatures  in  the  empira- 
cal  models  are  estimated  with  interior  ballistic  parameters  such  as 
charge  weight,  flame  temperature,  gun  diameter,  and  peak  pressure. ' 


2 

''Hypervelocity  Guns  and  the  Control  of  Gun  Erosion ,  "  Summary  Technical  ■ 
Report  of  Division  1 ,  NDRC,  Volume  1,  Washington ,  DC  1946. 

2 

J.A.  Lannon  and  J.R.  Ward,  "Workshop  Report  on  Wear-Reducing  Additives 
and  Their  Performance  in  Guns 3  "Proceedings  of  the  17th  JANNAF  Combustion 
Meeting3  CPIA  Publication  329,  November  1980. 

3 

A.J.  Bracuti3  L.  Bottei3  J.A.  Lannon 3  and  L.H.  Caveny3  "Evaluation  of 
Propellant  Erosivity  with  Vented  Chamber  Apparatus3  "  1980  JANNAF  Pro¬ 
pulsion  Meeting 3  CPIA  Publication 3  315,  March  1980. 

^ J.R.  Ward3  R.W.  Geene3  A.  Niiler3  A.  Rye3  and  B.B.  Gr oilman,  "Blowout 
Gun  Erosivity  Experiments  with  Double-Base 3  Triple-Base3  and  Nitramine 
Propellants 3  1980  JANNAF  Combustion  Meeting  CPIA  Publication3  315 3 
March  1980. 

5j. M.  Frankie  and  L.R.  Kruse 3  "A  Method  for  Estimating  the  Service  Life 

of  a  Gun  or  Howitzer,"  BRL  Memorandum  Report  No.  18523  June  1967. 

6  (AD  #818348)  . 

C.S.  Smith  and  J.S.  O'Brasky 3  "A  Procedure  for  Gun  Barrel  Life  Estima¬ 
tion,  "  Proceedings  of  the  Tri-Service  Gun  Tube  Wear  and  Erosion  Sympo¬ 
sium,  ADPA ,  Dover,  NJ,  March  1977. 

7 

J.  Comer,  Theori[  of  the  Interior  Ballistics^  of  Gunsj  John  Wiley  &  Sons, 
Inc.,  NY,  1950. 

Q 

L.W.  Nordheim,  H.  Soodak,  and  G.  Nordheim,  "Thermal  Effects  of  Propellant 
Gases  in  Erosion  Vents  and  Guns,"  NDRC  Armor  and  Ordnance  Report  No. 
A-262,  March  1944. 
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This  report  uses  Nordheim’s  scheme  to  compute  the  peak  surface 
temperatures  during  wear  tests  of  five  propellants.9 10  If  the  relative 
wear  of  the  propellants  in  the  blowout  gun  is  the  same  as  in  large  guns, 
then  a  plot  of  natural  logarithm  of  wear  vs  peak  temperature  should  be 
linear,  with  the  slope  identical  to  that  in  guns.  Wear  tests  ip  Reference 
9  incorporated  Niiler's  radioactive  technique^-9  to  correlate  mass  loss  and 
diameter  change. 

II.  INPUT  PARAMETERS  AND  COMPUTER  CODES 

Three  computer  codes  were  used  to  perform  the  calculations  dis¬ 
cussed  in  the  report.  The  first  program,  IB  3/70,  computes  interior 
ballistics  data  for  guns.  The  second  program,  NOZZLE,  computes  interior 
ballistic  data  for  a  blowout  gun  with  an  erosion  nozzle.  The  third  pro¬ 
gram,  NEWNSN,  computes  heat  transfer  using  files  of  IB  3/70  and  NOZZLE. 
Listings  for  each  program  are  contained  in  the  Appendix. 

The  assumed  heat  input,  as  described  fully  in  Reference  8,  was 
taken  to  be 


Q  =  JgXCppUAT  ,  (1) 

2 

where  q  =  heat  flux,  J/m  *s, 

X  =  friction  factor, 

C  =  specific  heat,  J/kg*K, 

3 

p  =  propellant  gas  density,  kg/m  3 
U  =  gas  velocity,  m/s,  and 

AT  =  temperature  difference  between  propellant  gas  and  wall,  K. 

Following  Nordheim's  method,  a  friction  factor  of  1/253  was  used 
for  the  17.3  mm  diameter  nozzle. 

Computations  were  done  for  the  five  propellants  evaluated  in  Re¬ 
ference  9.  The  compositions  are  listed  in  Tables  1  and  2. 


9  R.W.  Geene ,  J.R.  Ward,  T.L.  Brosseau ,  A.  Biller,  R.  Birhrtire,  J.J. 
Rocchio,  "Erosivity  of  a  Nitramine  Propellant, "  BRL  Technical  Report 
No.  02094,  August  1978.  (AD  #A060590) 

10S.E.  Caldwell  and  A.  Niiler,  " The  Measurement  of  Wear  from  Steel  Using 
the  Radioactive  hbCo ,"  BRL  Report  No.  1923,  September  1976.  (AD  #A030262) 
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TABLE  1.  COMPOSITIONS  AND  GRAIN  DIMENSIONS  OF  M30  AND  HFP  PROPELLANTS 


M30  HFP 

PPL-A-6372  PPL-A-6380 


Nitrocellulose  (12.6%N) 

28.0% 

29.3% 

Nitroglycerin 

22.5 

22.7 

Nitroguanidine 

47.7 

5.0 

RDX 

36.5 

Dioctylphthalate 

5.0 

Ethyl  Centralite 

1.5 

1.5 

Cryolite 

0.3 

Total  Volatiles  (Residual) 

0.2 

0.3 

Grain  Length,  mm 

7.78 

10.58 

Grain  Diameter,  mm 

1.59 

2.37 

Grain  Perf.  Diameter,  mm 

0.46 

0.77 

Grain  Web,  mm 

0.56 

0.80 

Grain  Geometry 

SP 

SP 

TABLE  2.  COMPOSITIONS  AND 

GRAIN  DIMENSIONS 

OF  M5 ,  M8 ,  AND  Ml 

PROPELLANTS 

M5_ 

M8 

Ml 

Nitrocellulose  (13.25%N) 

81.95% 

52.15% 

85.00% 

Nitroglycerin 

15.00 

43.00 

Ethyl  Centralite 

0.60 

0.60 

Barium  Nitrate 

1.40 

Potassium  Nitrate 

0.75 

1.25 

Diethylphtalate 

3.00 

Dinitrotoluene 

10.00 

Dibutylphthalate 

5.00 

Diphenylamine,  Added 

1.00 

Ethyl  Alcohol,  Residual 

2.30 

0.40 

0.75 

Water,  Residual 

0.70 

0.50 

Graphite 

0.30 

Grain  Length,  mm 

10.58 

25.4 

8.26 

Grain  Diameter,  mm 

3.92 

12.7 

3.68 

Grain  Perf.  Diameter,  mm 

0.41 

0.37 

Grain  Web,  mm 

0.69 

0.56 

0.64 

Grain  Geometry 

7  Perf 

Strip 

7  Perf 
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Propellant  parameters  required  for  the  computations  are  listed  in 
Table  3. 

The  term,  B/W,  is  used  to  relate  burning  rate,  r,  with  chamber  pressure,  P, 
and  propellant  web,  W,  by 

r  =  (P  +  44.8) (|)  •  W  .  (2) 

Propellant  webs  are  listed  in  Tables  1  and  2. 

3 

For  the  blowout  gun,  chamber  volume  and  nozzle  area  were  315  cm 
and  2.36  cm2,  respectively. 


III.  RESULTS  AND  DISCUSSION 

The  experimental  results  from  Reference  9  are  listed  in  Table  4. 

Table  5  gives  calculated  time  to  maximum  pressure,  peak  nozzle  sur¬ 
face  temperature,  heat  input,  and  time  from  nozzle  rupture  to  peak  nozzle 
surface  temperature. 

A  linear  least  squares  fit  tested  the  exponential  dependence  of  wear 
with  peak  surface  temperature. 

£n  (w)  =  a  +  bT  ,  (3) 


where  w  =  wear,  y/shotj 

T  =  peak  bore  surface  temperature,  K,  and 
a,b  =  constants. 

The  best  fit  values  of  a  and  b  were  -9.7  ±  2  and  0.0094  ±  0.002, 
respectively,  with  error  representing  one  standard  deviation.  Table  6 
compares  experimental  wear  and  wear  computed  with  Eq.  (3).  Eq.  (3) 
clearly  underestimates  wear  for  Ml  with  a  peak  temperature  of  773  K. 
Frankle-Kruse  also  noted  their  expression  was  invalid  below  900  K; 
Smith-0' Brasky  report  a  similar  threshold  of  750  K. 

A  check  on  how  wear  varies  with  peak  bore  surface  temperatures  in 
guns  can  be  made  with  bore  surface  temperature  computations  by  Vassallo 
and  coworkers  at  Cal span  Corporation! 1  in  the  M185  cannon.  The  Cal span 
workers  compute  peak  bore  surface  temperatures  from  measured  total  heat 
input. 12  Recent  work  in  a  shock-tube  gun  with  different  metals  showed 
temperature  calculations  starting  with  measured  total  heat  input  could 
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TABLE  3.  PROPELLANT  PARAMETERS* 
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TABLE  4.  EXPERIMENTAL  DATA  FROM  BLOWOUT  GUN  EXPERIMENTS 
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TABLE  5.  COMPUTED  HEAT  TRANSFER  RESULTS 
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TABLE  6 


COMPARISON  BETWEEN  EXPERIMENTAL  WEAR  AND  WEAR  COMPUTED 
WITH  EQUATION  (3) 


Propellant  Peak  Temp. ,  K  Wear  exp't,  y/Shot  Wear  Cal’d,  n/Shot 


Ml 

Ml 


773 

875 


0.2 

0.1 


0.0005 

0.2 


M30 

966 

0.5 

M30 

1,130 

0.6 

M30 

1,228 

3.8 

0.5 

2.5 

6.2 


HFP 

968 

0.5 

HFP 

1,150 

1.1 

HFP 

1,235 

6.8 

0.5 

3.0 

6.6 


M5 

982 

0,8 

M5 

1,155 

4.1 

M5 

1,250 

18.5 

0.6 

3.1 

8.3 


M8 

1,107 

M8 

1,210 

MS 

1,407 

2.8 

9.6 

48.6 


2.0 

5.0 

32.1 
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predict  correctly  whether  the  melting  point  of  the  metal  was  exceeded. 
Starting  conditions  in  the  shock-tube  gun  were  varied  to  produce  tempera¬ 
ture  bracketing  the  metals'  melting  temperatures. 

Table  7  shows  bore  surface  temperatures  computed  for  three  155mm 
propelling  charges  which  have  no  wear-reducing  additive.  Wear  data  is 
available  for  the  M119  and  the  XM201E2  charges, 14 >15  Wear  sensor  data  in 
the  M185  cannon  showed  the  XM208  minus  its  additiye  was  five  times  as  ero¬ 
sive  as  the  XM201E2.  Wear  of  the  XM208  was  presumed  to  be  five  times 
greater  than  the  measured  wear  of  the  XM201E2  charge.  A  least-squares 
fit  of  the  three  points  in  Table  8  to  Eq.  (3)  produced  values  of  a  and  b 
of  -7.4  and  0.0076K-!,  respectively. 

TABLE  7.  BORE  SURFACE  TEMPERATURES  AND  WEAR  FOR  155mm  CHARGES  WITHOUT 
ADDITIVE 


Charge  Peak  Temp,  K  Wear,  p/shot  Ref. 


M119 

961 

0,9 

15 

XM201E2 

1,100 

2.6 

16 

XM208  (no  additive) 

1,306 

12.5 

12 

Figure  1  plots  wear  vs  peak  bore  temperature  for  the  nozzle  and  the 
155mm  gun.  Figure  1  suggests  relative  wear  among  propellants  will  be 
overestimated  compared  with  wear  in  guns  by  £en  percent.  Nonetheless, 
the  agreement  seems  close  enough  to  suggest  relative  propellant  erosivity 
can  be  inferred  in  nozzles. 
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Figure  1.  Wear  vs_  Peak  Bore  Surface  Temperature 
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Some  other  observations  from  the  calculations  of  bore  surface  tem¬ 
peratures  can  be  made  from  comparisons  of  bore  surface  temperature  vs. 
time  for  the  37mm  blowout  gun  and  a  3"/70  Naval  gun  depicted  in  Figures 
2-4.  Table  8  summarizes  pertinent  heat  transfer  data  for  the  three 
figures . 


TABLE  8.  CALCULATED  HEAT  INPUT  AND  SURFACE  TEMPERATURES  FOR  3"/70 
AND  37mm  BLOWOUT  GUNS 


Gun 

Nozzle 
Diameter , mm 

1/A 

Propellant 

Charge 

Mass,g 

Flame 

Temp,K 

Surface 
Temp , K 

Heat  In-^ 
put, J/mm 

37mm 

17.3 

253 

HFP 

89.8 

3,255 

1,235 

0.78 

3"/70 

- 

280 

Picrite 

3991. 

2,065 

1,224 

1.13 

37mm 

13.3 

239 

Ml 

159 

2,417 

1,254 

1.21 

The  first  point  of  interest  is  the  short  heating  time  of  the  blow¬ 
out  gun  in  Figure  2  vs.  the  3"/70  Naval  gun.  As  a  result  of  the  longer 
heating  time,  the  3"/70  gun  reaches  the  same  peak  bore  surface  tempera¬ 
ture  with  a  much  cooler  propellant.  The  effect  of  the  longer  heating 
time  is  also  reflected  in  the  larger  total  heat  input  for  the  gun.  The 
figures  also  point  out  why  propellants  with  flame  temperatures  near 
3800K  would  be  needed  to  get  the  wear/shot  characteristic  of  tank  guns 
with  the  17.3  mm  nozzle.  Figure  4  shows  that  Ml  propellant  produces  an 
equivalent  peak  temperature  as  HFP  if  charge  mass  is  increased  and  nozzle 
diameter  is  decreased.  A  new  chamber  would  be  required  in  the  37mm  gun 
to  be  able  to  load  160g  of  propellant. 


IV.  CONCLUSIONS 

1.  Bore  surface  temperature  vs.  time  were  computed  for  17.3  mm  diameter 
nozzles  in  the  BRL  37mm  blowout  gun  in  order  to  test  dependence  of  wear 
vs.  peak  bore  surface  temperature. 

2.  The  natural  logarithm  of  wear  from  the  nozzle  was  shown  to  be  linearly 
dependent  on  peak  bore  surface  temperature  which  agrees  with  empirical 
wear  models  and  with  calculations  of  peak  bore  surface  temperatures  by 
Calspan  Corp.  for  the  155mm  M185  cannon.  It  was  shown  the  slope  is  10% 
steeper  in  the  nozzle  suggesting  relative  propellant  erosivity  with  the 
17.3mm  diameter  nozzle  may  be  exaggerated  slightly  compared  to  large 
caliber  guns. 
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Figure  2.  Temperature  vs.  Time  In  37mm  Blowout  Gun  (17.3mm 
Diameter  NozzTe-HFP  Propellant). 
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TEMPERATURE  (DEG.  K ) 
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Figure  3.  Temperature  vs.  Time  In  3"/70  Gun  With  Plcrite  Propellant. 
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TEMPERATURE  (DEG.  K) 
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Figure  4.  Temperature  vs.  Time  in  37mm  Blowout  Gun  (13.3mm 
Diameter  NozITe-Ml  Propellant). 
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APPENPIX 


COMPUTER  CODRS 


COMPUTER  CODES 


The  first  program,  "IB3/70",  computes  interior  ballistics  data 
for  guns. 

The  second  program,  "NOZZLE",  computes  interior  ballistics  data 
for  a  blowout  gun  with  an  erosion  nozzle. 

The  third  program,  "newnsn"  computes  heat  transfer  and  conduction 
data.  The  program  "newnsn"  uses  information  stored  on  files  "DATA" 
and  "DATA1"  by  "IB3/70"  or  "NOZZLE". 
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"IB3/70" 


Input  data  begins  in  line  2860.  Data  to  be  stored  in  each  line 


is  as  follows: 

Line 

Data 

2860 

Projectile  weight  (lbs) 

2870 

Projectile  travel  (inches),  chamber  volume  (cu.in.) 

2880 

Pidduek-Kent  constant,  integration  step  size  (sec) 

2890 

Estimated  projectile  velocity  (ft/sec),  bore  diameter (in) , 
bore  area  (sq.in.) 

2900 

Shot  start  pressure  (psi) 

2910 

Igniter  weight  (pounds) 

2920 

Igniter  impetus  (in  lbf/lbm),  ratio  of  specific  heats, 
flame  temperature  (°K),  covolume  (cu.in. /lb) 

2930 

Propellant  weight  (lbs) 

2940 

Propellant  force  (  in  lbf/lbm) ,  ratio  of  specific  heats, 
flame  temperature  (°K),  covolume  (cu.in. /lb),  density (lb/cu. in 

2950 

Propellant  grain  diameter  (in),  length  (in),  perforation  dia¬ 
meter  (in) ,  and  number  of  perfs 

2960 

Propellant  identification  (an  alphanumeric  string) 

2970 

PI,  B (where  the  burning  rate  (in/sec)  is  B* (space  mean 
pressure  in  psi  +  PI)),  number  of  integration  steps  per 
printout,  expected  peak  chamber  pressure  (psi),  expected 
peak  acceleration  (g's). 

2980 

Number  of  "resistive  pressures'1  listed,  followed  by  the 
appropriate  number  of  pairs  of  travel  distance (in)  and 
resistive  pressure  (lb/sq.in.) 

2990-3020 

These  lines  are  used  by  the  program  to  aid  in  calculating 
the  fraction  of  the  7  perf  grain  burnt  after  splintering. 

3025 

Use  this  line  to  provide  data  if  called  for  by  the  input 
statements 

Data  requested  by  input  statements: 

"Number  of  parameters  to  be  varied":  Enter  a  number  (0  to  5). 

"Number  of  proj  wt,  cham  vol,  pmax,  travel  distance,  charge  wt":  The 
number  of  non-zero  enteries  should  be  as  indicated  in  the  immediately 
previous  input.  The  values  for  the  parameters  indicated  should  have  been 
listed  in  line  3025  of  the  program. 

"Option  No  2" 

Enter  a  number  (0  to  7).  The  program  will  then  perform  as  indicated 
below: 
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NUMBER  ENTERED 


ACTION 


0  Compute  IB  trajector  for  input  data  (i.e.,  data 

in  lines  2860-2980) 

1  Change  B  to  match  indicated  peak  pressure 

2  Change  charge  weight  to  match  indicated  pressure 

3  Change  B  to  match  indicated  acceleration 

4  Change  charge  weight  to  match  indicated  acceleration 

5  Change  charge  weight  to  match  indicated  velocity 

6  Print  out  a  blurb  of  the  above 

7  Same  as  0,  but  program  writes  data  on  files  DATA 
DATA1  for  future  use. 

Note  that  if  the  first  input  requested  was  not  zero  multiple  runs  will  be 
performed.  This  is  not  useful  if  "Option  number  2"  is  7. 
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10. 

2  0  0 
30 
40 
50 
60 
?0 
SO 
90 
100 
110 
120 
130 
140 
150 
160 
17G 
ISO 
190 
2  0  0 
210 
220 
236 
246 
250 
2  6  0 
270 
2  S  0 

2  9  0 

3  O  0 
310 
3  2  O 
330 
3  4  0 
350 

3  6  6 
370 
380 
390 

4  00 
418 
420 
430 
440 
450 
460 
470 
480 
496 

5  0  0 
510 
520 
536 
540 

559 

560 
570 

580 

581 
598 
600 
6 1  O 
620 
630 
64  0 
650 


OPTION  BRSE  1 
VtRLHP 
DIM  C < 5 0 > 

DIM  F^C 103 

HIM  T  <  2 ,  1 0  ,  F  <  2 ,  1  2  )  ,  R  <,  6 , 3  ) 

Eheat *=0 
LET  Z6=0 

D  E  F  F  H  fl  <  Y  9  >  =  I H  T  <  Y  9  *  1 0  6  0  >  / 1  O  0  0 
R  E  R  D  L  1  ,  G  n  L>  3  ,  G  5 ,  L,  6  ,  G  7  ,  0  3  j  G  4  ,  G  9 
LET  E7=G?*12 
READ  11,12,13,14,15 

RE AH  P 1 , P2 , P 3 , P4 , P5, P6, P7, P8, PS, Py 
READ  PiT 

READ  H ,  E ,  M 2 ,  M  3 ,  A  9 
IF  P0< >7  THEN  190 
LET  M  1  =  <  P  7  - 3 * P 9  )  /■  4 

LET  W2=<  <P7-W1  >  A2-?*<jP9+Wi  ) 2  > / < P ? *  2  -  7  *  P 9 -•  2 > 

GOTO  200 
LET  W 1 ~ 1 0 0 0 
READ  K 

FOR  1=1  TO  K 
READ  TO,  I  >  ,  T  .<  2 ,  I  > 

NEXT  I 

FOR  1=1  TO  12 
READ  F(1JI),F(2,() 

NEXT  I 

PRINT  " nUMBER  OF  PARAMETERS  TO  PE  VARIED" 

INPUT  D 

IF  D = 0  THEN  510  * 

PRINT  "NO  OF  -PRO J  NT, CHAM  VOL , PMAX, TRAVEL  HIST, CHARGE  WTM 

INPUT  D 1 , D2 , D3 , D4 , D5 

IF  Dl»©  THEN  366 

FOR  N 1 = 1  TO  D1 

READ  L<N1>. 

NEXT  Nl 

IF  D 2  =  0  THEN  406 
FOR  N 2=1  TO  D2 
READ  MCN2> 

NEXT  N2  ; 

IF  D3=0  THEN  470 
FOR  N 3  = 1  TO  D3 
READ  0 < N 3 > 

NEXT  N3 

FOR  N 4  = 1  TO  D4 
READ  W <  N4  ) 

NEXT  N4 

IF  D 5 = 0  THEN  510 
FOR  N 5  =  Q  TO  D5 
READ  YCN5> 

NEXT  N5 

PRINT  "OPTION  NO  2  (0<^X>»5>  ;  INPUT  6  fOR  EXLANAT I ON " 

INPUT  Cl 

IF  Cl <6  THEN  760 
IF  Cl =6  THEN  620 
C1=0 

CREATE  " DAT  ft " ,  126 
CREATE  " DATA 1 “ , 3 
ASSIGN  #1  TO  "DATA" 

BUFFER  #1 

ASSIGN  #2  TO  " D ft T ft 1 " 

Eh eat = 1 
GOTO  930 
PRINT 
PRINT 
PRINT 

PRINT  "TO  MATCH";  TAB  <20  );  "BY  ADJUSTING'1 ;  TfiB<4  5  11  INPUT" 
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660 

PRINT 

670 

PRINT  “NOTHING" ; TREK  4 

680 

PRINT  "  PMAX" 

;  TAB*: 20) ; 

698 

PRINT  "PM  FIX" 

; TflB<28> ; 

700 

PRINT  11  ftCCEL 

" ;  T  ft  B  <  2  0  > 

716 

PRINT  "  ftCCEL 

"  ; TAB <20 ) 

720 

PRINT  "VELOC 

I  TV 11 ;  TRB< 

738 

PRINT 

740 

PRINT 

750 

INPUT  Cl 

760 

IF  B = O  THEN 

930 

770 

IF  Ii  1=8  THEN 

800 

f'  y  0 

FOR  N 1 = 1  TO 

D1 

790 

LET,  6 1  =  L  <  N  1  ) 

80G 

IF  B 2  —  0  THEN 

830 

8  1  0 

FOR  N2= 1  TO 

D2 

" CHARGE 


WT,,;TflB<43>| 

T  Ci  V.  ( *  A  O  s,  •  «•  O  I' 


20)  ;  •' CHARGE  NWT"  ;  T  ft  E  <  4 3 ) 


820  LET  G 3  =  M < N 2 > 
830  IF  H3=G  THEM 


;i  6  0 


848 

FOR 

H3= 1  TO  D3 

850 

LET 

M  3  Q  (  H  3  > 

860 

IF  D 4  =  0  THEN  890 

878 

FOR 

N4= 1  TO  B4 

8  y  u 

LET 

G2  =  H  <  N4 ) 

890 

IF  D5=0  THEN  930 

900 

FOR 

N5=i  TO  B5 

91G 

LET 

P  1  =  Y  <  N  5  ) 

9  2  0 

NEXT  N5 

930 

LET 

J  1  =  I 1*12/ < 13-1 > 

940 

LET 

J  2  =  J  1  /  I  4 

950 

LET 

J3= 11*1 2/ I 4 

9  G  O 

LET 

J  5  = I 1*15 

978 

LET 

Q 1  =  P 1  *  P 2 / < P  3-1 ) 

980 

LET 

Q  2 = Q 1  /  P  4 

9  9  0 

LET 

Q  3  =  P 1  *  P  2  /  F*  4 

1  G00 

LET 

Q4  =  P  1  /F*6 

1010 

LET 

Q  5  =  P 1  *  P  5 

1020 

LET 

Q6=I 1+P1 

1030 

Q7  = 

P  8  *  <  P  7  •  2-P0*P  9  2  ) 

1040 

LET 

Q  8 = G  6  /  2 

1  050 

LET 

0 1  =  ■  y  8  *  G  8  1 . 5  *  <  G  2 + G  3  /  G  4 )  * 

1  060 

LET 

0 1  =0 1  /  <  1  + .  S*G8A2 .  175/P  1 A . 

1070 

LET 

H 1 =0 1  +  C  G 1 +Q6/G5 ) /?72 

1  080 

LET 

0  2=1  /  <  P  3  - 1  ;• 

1090 

LET 

H2= 1 +Q6/<G 1 *G5 ) 

1100 

LET 

03= 1 / <  < 2*02  +  3 ) /G5+2* <  02+ 1 

1110 

LET 

H3=<  1-03)  M-02-1  ) 

1  120 

LET 

H  4  =  G  4  *  3  8  6  /  G 1 

1130 

IF 

M2>= 1 50  THEN  1160 

<  p  4  -  3  0  0  > 

8  3  7  5  >  /  <  1  4  4  *  G  7  *  G  7  ) 


+  1  >*61/06) 


1140  PRINT  "TIME  BIST 

1150  PRINT  "FR.  RC CEL 

1160  LET  T  = 1 0 0 0 
1170  LET  T  2  =  T  <  2 ,  1) 

1180  LET  H  =  0 
1190  LET  M  1  =  1 

1200  LET  R  (  1  ,  1  )  =  R  <  1 , 2  )  =  R  <  1  , 
1210  LET  X 1  =  R  <  1 , 2 > 

1220  LET  V 1 =R  < 1 , 3) 

1230  FOR  1=1  TO  3 
1240  LET  R<2,  I  )  =  R < 1,1) 

1250  NEXT  I 

1260  LET  J  =  2 

1278  GO SUB  2280 

123G  LET  S-S+Q8 

1290  LET  J  =  3 

1300  FOR  1=1  TO  3 

1310  LET  R<2, I)=R<2, I )*G6 


PRES 

BRSE" 


VEL 


TEMP 


>=B0=Z"T2=S=B3=B4=B5=Z5»Z6*0 
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1320  LET  R  •;  3  ,  I  )  -  R  <1,1  >  +  R  C  2 ,  I  >  **'  2 

1336  NEXT  I 

1  3  4  0  G  0  S  U  B  2  2  8  8 

1350  LET  J =4 

1360  FOR  1=1  TO  3 

1370  LET  R ( 3 ,  I >  =  R  <  3 , I >  * G 6 

1380  LET  R  <  4  ,  I  >  =  R  <  1  f  I  >  +  R < 3 ,  I  )  /*2 

1390  NEXT  I 

1  4  0  0  G  0  S  U  E  2 2  8  O 

1410  LET  J=5 

1420  LET  S=S+Q8 

1430  FOR  1=1  TO  3 

1440  LET  R<4, I )=R<4, I )*G6 

1450  LET  R <  5  .  I  >  =  R  <  1  ,  I  >  +  R  <  4  ,  I > 

1460  NEXT  I 

14  7 0  GOSU B  2280 

I486  FOR  1=1  TO  3 

1490  LET  R<5, I >=R<5, I >*G6 

1500  LET  R ( 1 , I >  =R < 1 , I)  +  <R<2, I >  +  R < 3 , I)*2+2*R<4? I>+R<S, I >>^6 
1510  LET  R<6, I >=R< 1 ,  I  > 

1520  NEXT  I 

1530  LET  ,T=6  s 

1540  IF  Eh  eat  THEN  PRINT  #1 ;  V,  X,  Pl*ZXG3+G4*)<-Pi*U~Z>/P6>  ,  J 
1550  IF  M2*  I  NT  <  M 1  /M2  >  <  >P1 1  THEN  1590 

1 560  PR  1  NT  S* 1 G00 ; TflB < 8 > ; FNR (R(l,2)); TAB < 1 8  > ; I NT  < BQ > ; TAB < 27 > 5 
157G  PRINT  I  NT  <  R  <1 ?  3  > / 1 2  > ; TAB  < 34 ) ; I NT < T  > ; TAB <45> ; FNft  < Z > 5 
1580  PRINT  TAB < 54); I  NT < A 1 /386 . 09 > ; TAB < 63 > ; IHT<B2> 

159U  IF  R < 1 , 2 ) >  =  G 2 - 2 5  THEN  1610 

16G0  GOTO  1638 

1610  H  =  N  + 1 

1626  LET  C < N >=B2 

1630  IF  R < 1 , 2) >=G2  THEN  1730 

1640  LET  M 1  =  M 1  + 1 

165G  LET  V  =  R< 1,3) 

1668  IF  B 0 <  = B 3  THEN  1210 
1670  LET  A  5  =  A 1  / 3 8 6 . 0 9 
1630  LET  B3=B8 

1690  LET  B 8  =  B 2  j 

1700  LET  B 4  =  8 
1718  LET  E 5  =  R < 1 , 2 > 

1720  GOTO  121© 

1730  LET  X  =  <  G  2  - X 1 ) / < R<1,2)~X1> 

1748  LET  V  =  V 1 +  X * < R(l, 3 > - V 1 > 

1750  LET  3 1 =  S - Q 8 

1760  LET  Z9=C<N-i>+X*<C<N>-C<N-l>> 

1770  IF  C 1=0  THEN  I960 
1  7  8 0  B r  an c  h = C 1 

1798  IF  < Bnanc h >= 1 )  AND  <Branch<6)  THEN  ON  Branch  GOTO  1800, 1830, 1868, 1898, 1920 
1800  IF  ABSCB3-M3X200  THEN  I960 
1810  LET  B= B s  < 1 - . 85*  <  M 3 - B 3  }  /  <.  M3 + B 3  )  ) 

1320  GOTO  1940 

1830  IF  ABS<  B3-M3X200  THEN  I960 
1840  LET  P 1 =P 1 / <  1  ~ . 9*  < M3-B3 > / < M3+B3 > > 

1850  GOTO  1940 

I860  IF  A  B  S  <  A  9  -  A  5  >  < 10 0  THEN  I960 
1870  LET  B  =  E s ( 1-1  . 1  * < A9- A5 > ✓ < A9  +  A5 > ) 

1880  GOTO  1940 

1890  IF  A  £  S  <  A  9  -  A  5  )  < 1 0 0  THEN  I960 
I960  LET  P  1  =P  1  ✓  1-.  9  *  <  H  9  -  A  5  >  /  <  A  9  +  A  5  >  ) 

1910  GOTO  194G 

1920  IF  HB3(G7“V/Q2)< 10  THEN  I960 

1938  LET  P1=P1X1-1 .5*<G7~V/pi2>XG?+V<'12>> 

1940  LET  E7=V 
1950  GOTO  930 
I960  PRINT 

1970  IF  Eheat=0  THEN  2000  29 


1980 
2000 
2010 
2O20 
203O 
2040 
2O50 
2060 
207O 
2030 
2090 
2100 
2110 
2120 
2130 
2  1  40 
2150 
2160 
2170 
2100 
2  1  90 
220O 
221Q 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
2340 
2358 
2360 
2370 
2380 
2390 
2400 
2410 
2426 
2430 
2448 
2458 
2460 
2470 
2480 
24  90 
2500 
2510 
2526 
2530 
2540 
2550 
2560 
2570 
2580 
2590 
2600 
2610 
262G 
2630 
2648 


CALL  Ej  C#1  ,  #  2  ,  P  1  .*•'  <  G  3  +  G  2  *  G  4  >  ,  <  V  >  ,  S,  G6,  T  ,  P5f  G2  +  G3/G4,  G8,  G3.'G4  > 

IF  M2 >=500  THEN  2190 

PRINT  " PRO J . NT .  = " ; G 1 ; "LBS . " 

PR  I  NT  “  I  HI  T  .  CHAM .  VOL .  =  "  ;  G3 ;  "  CU  -  I N .  " 

PRINT  " CHG  NT  ; PI ; “LBS. " 

PRINT  “TRAVEL  BIST-  = “ ; G2 5 “ I N . " 

PRINT  "MUZZLE  VELOCITY  =  "  ;  I  NT<  V/ 12);  “F*'S  “ 

PRINT  "MAN.  CHRM.  PrES.  =  “  ;  I  NT  <  B3  >  ;  "  PS  I  flT“ ;  JES4;  "SEC-  OR"  ;  B5; 11 JN,  " 
PRINT  "MAX.  BASE  PRES.  * " ; I NT < B8> ; “PS I" 

PRINT  "MUZZLE  PRESSURE  =  “ ; I  NT < Z9 > ; “ PS I  AT  M ;  SI  J “ SEC" 

PRINT  "MAX.  ACCEL  = “ ; I  NT  < A5  > ; " GS " 

IF  Z<1  THEN  2130 

PRINT  "FRACTION  BURNT  = “ ; 2: “AT" ; 25; " IN  FROM  MUZZLE" 

GOT 0-  2140 

PRINT  "FRACTION  BURNT  =" ;2 

PRINT  "BETA" ; B 

PRINT 

PRINT  G9 s I  1 , T < 1 , 1 > , T < 2 , 1 > , T < 1 , 2) , T < 2, 2> 

PRINT 

GOTO  2740 

IF  E9= 1  THEN  2240 

PRINT  " PRO J  NT  BIST  CHAM  VOL  CM  PMAKP 

PRINT  "IV  ACCEL  BETA" 

LET  E9=  1 
PRINT 

PR  I  NT  TAB  C 1 > ; G 1 ; TAB  < 8 > ; G2 ; TAB < 1 5 >  J G3 ; TAB < 24  >  ;  P  1 ; TAB < 3£  > : 

PR  I  NT  I  NT  < B3 > ; T ; B <48 > ; I  NT  < V/ 1 2 > ; TAB  <48 > ; I NT <  A5  > ; TAB <  56 > ; B 

PRINT 

GOTO  2740 

LET  U  =  R < J ,  1 >  • 

LET  X  =  R  <  J  ,  2 ') 

LET  V  =  R  <  J , 3 ) 

IF  2>= 1  THEN  2350 
IF  U >  =  M 1  THEN  2400 

LET  Z  =  1  -  <  P8-U  )  *  <  <  P7-U  >  A-2 -PO*  <  F9  +  U  )  -*  2  >  /Q7 
GOTO  2470 
LET  2=1 

IF  26=1  THEN  2470  j 

LET  2 5 = G 2 - X 

LET  26=1 

GOTO  2470 

FOR  K = 1  TO  12 

IF  FC1,K>*W1>U  THEN  2440 

NEXT  K 

GOTO  2350 

LET  M  3  = < U/W i -F  C 1 , K- 1 >  > ✓ C F  <1 , K  >  ~F  <1 , K- 1 > > 

LET  M 4  =  < P 8 - U > / P 8 

LET  2=1-U4*<FC2?  K-1>+CF<2, K>-FC2P K-l > >*W3>*W2 
LET  T 1 =0 

IF  X  =  0  THEN  2560 

FOR  K=2  TO  10 

IF  X <  =  T < 1 , K >  THEN  2530 

LET  T 1 =T 1  +  < T  < 1 , K ) -T <l,K-l>)*<TC2,IO  +  T<2, K- 1 >  > ,'2 
NEXT  K 

LET  T2  =  T  <  2 , K- 1 >  +  <  T  <  2 , K > -T  < 2 *  K- 1 > ) *  <  T  < 1 , K > -X ) / C  T  < 1 , K > -T  < 1 , K~ 1 >  > 
LET  T 1 =T 1  +  <  T2  +  T  <  2 , K- 1 >  > * < X-T C 1 ,  K-l > > s2 
LET  T1=T1*G4 

LET  T=CJ1+Q1#Z-H1* V * V - T 1 >/( J2+Q2*Z> 

LET  V0  =  G3  +  G4*X-Q4#  < i  -2  > - J5-Q5*Z 
IF  V 0 < 0  THEN  2680 
GOTO  2620 

PRINT  G3,G4,Q4, J5,Q5.Z 
GOTO  3030 

LET  B 1 =  T  *  <  Q  3  *  Z + J  3 ) ✓ V 0 
LET  B 2  =  B 1 / H 2 

LET  B0=B2*H3  30 


t 


I 


* 


2656 
266  U 
26?G 
2680 
2690 
2708 
2710 
2720 

2739 

2740 
2750 
2760 
2770 
2788 
2790 
2800 
2810 
2820 
2830 
2840 
2850 
2868 
2370 
2880 

2899 

2900 
2910 
2920 
2930 
2940 
2950 
2960 
2970 
2  9  8  0 
2998 

3009 

3010 
3020 
303O 
3998 
4000 
4005 
4020 
406O 
4070 
4120 
4130 
4140 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
422G 
4230 
4240 
4250 
426u 
4  270 

4280 

4281 
4290 

4300 

4301 


IF  X<>8  THEN  2670 
IF  B1CG9  THEN  2S90 
LEI  A1=H4*<B2-T2> 

GOTO  2700 
LET  A  1=8 

LET  R< J, 1 >=2*B*<B1+R> 

LET  R < J , 2)= V 
LET  R < J , 3  >  =  fi  1 
RETURN 

IF  H=0  THEN  2850 
IF  D5=0  THEN  2770 
NEXT  N5 

IF  H 4 <  =  0  THEN  2790 
NEXT  N4 

IF  D 3  =  0  THEN  2819 
NEXT  N3 

IF  D2=e  THEN  2830 
NEXT  N2 

IF  B 1  =  0  THEN  2850 

NEXT  HI 

PAUSE 

BATfi  15 

BATA  185.84,360 

BATA  3,. 0O005 

BATH  3400, 3,7.21 

DATA  2500 

BATA  .02 

BATA  1 . 152E+06, 1 . 25,2000,30 
BATA  8.8 

BATA  3445O00, t . 269,2065,32.81 , . 056- 
BAT  A  . 1425, . 35, .0132,7 
DATA  11 EX6586" 

DATA  6500, . 0 0 0 0 56673, 200,52380,0 
DATA  2,0,1808,508,1 888 

BATA  1,1,1 .006, 8.  967,  !  .  024, 0. 874,  1 . 05.5, 8.735 
BATA  1.101,8, 569, 1 . 162,0.395, 1 . 24,8.233 
BATA  1 . 292,0. 191,1. 347,8. 104,1. 489, 0,049 
BATA  1 . 48,0.014, 1 . 557,0 
STOP 

SUB  E j  <  # 1 , #2 , Rho i , V i , T i , Bt , T erop i , Ada, Xm . B , C 1 > 
OPTION  BASE  1 
OVERLAP 
E  o  a=  E  o  b  9 
X  t  =Xm 

Kount =  T i /lit 
V i = V i /Xm 
Rh  =  1  /'Rho  i 
T  p  =  0 
Gm 1 = . 25 
Cm=. 34 

I NPUT  " C v , Hda, Gm 1 " , Cv , Ada, Gm 1 

A 110: K o u n t = K o u n t  + 1 

T  p  =  T  +  D  t 

T=T i +Tp 

Bem=l +Vi *Tp 

V  =  V  i  /  li  e  m 

Rho = Rho  i  .'Bern 

Ternp=T emp  i  *<  <Rh-fldaV<  1  /Rho- Ada)  ) A Q m  1 
V=V*Xm 

PRINT  #1 ; V,Xt, Rho, Temp 
IF  Kount  MOB  408  THEN  Ally 
IF  Eoa  THEN  A40 
IF  V*Rho>=l©0  THEN  A110 
Eoa=  1 
D  t  =  4  *  D  t 

PR IHT  #1 ; 1E51 ,0,0,0  31 


CO  CO  CO 


4305  GOTO  ft  110 

4310  ft40:IF  Rho>=. 00050  THEN  ftl  10 
4315  PRINT  #1 ;  1E51 , 0, 0, 0 
4320  Km=Xt 
4330  Dt=Dt/'4 

40  PRINT  #2;  Rhoi  ,  Vi  ,Ti  ,  fif,  Tempi  ,  ftda,  Xm ,  I',  C 1 
60  SUBEND 
70  END 


t 
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J 


"NOZZLE" 


Input  data 
of  lines  20  and 

Parameter 

F0 

Beta,  PI 

C 

Rhop 

Ada 

TO 

Gamma 

Force 

C0,C2,C2P 

VO 

Ast 

Pstart 

G6 


is  lines  110  and  120  corresponds  to  the  READ  statements 
30.  Input  dat  is  as  follows: 

Data 

Initial  web  fraction  remainging.  (0<F0<1) 

If  F0  is  too  small  the  initial  pressure  will 
exceed  the  intended  nozzle  start  pressure. 

The  assumed  burning  rate,  in  inches  per  second, 
is  Beta  (P+P1)/W  where  P  is  the  chamber  pressure 
in  psi  and  W  is  the  web  size  (inches) 

Charge  weight  (pounds) 

Propellant  density  (lb/cu.in.) 

Propellant  covolume  (cu. in/lb.) 

Adiabatic  flame  temperature  (°K) 

Ratio  of  specific  heats 

Impetus  (in  lbf/lbm) 

These  are  the  kp,  k2,  and  -k2  of  NDRC  A142  of 
5  Feb  1943;  they  are  used  to  help  calculate  the 
amount  of  propellant  that  has  been  burnt. 

Chamber  volume  (cu.in.) 

Area  of  the  erosion  nozzle  (sq.in.) 

Nozzle  start  pressure  (psi) 

Integration  step  size  (sec) 
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10  OVERLAP 

2  0  R  E  A  D  F  0 ,  B  e  ta,Pl,C,  R  h  op,  Ada,  T  0 ,  G  am  rn  a,Fo  ►'  ce,  CO.  C  2 ,  C  2  p 

SO  READ  V8 , Hst , Pst  art , G6 
40  S i =S=0 

50  2eta=(2/(Gamma+l  )  (  1  /(Gamma- 1  )  ) 

60  Ep-< 1 - 2  *  Z  a  t a  )  /  G  am  m  a 

70  P  h  e  1  p  =  ( 2  /  (  G am m a+ 1 > ) A < G am m a/ < G am m a- 1 >  > 

80  Q 8=66/2 
90  Break  =0 

1  OS  I  MAGE  DD . BOLD , DD . DDL , DDDD . DDUD , DDDDDD , DDDDDDDDD , PHD. BDBBB , 8D 
j 1 0  IlRT fl  ®  •  500 ,  .  005200 , 6500 ,  .  300 ,  .  05? ,  32 .81, 2065 ,  1 . 269 , 3445000 ,  .  8502 

120  DATA  19. 2, . 098, 41000, . 000050 

130  R  <  1  1  )  =  F  0 

140  R(l,2>=0 

150  Kount=0 

160  CREATE  “DATA", 126 

170  CREATE  " DA T A 1  ",  3 

180  ASSIGN  #1  TO  "DATA" 

190  ASSIGN  #2  TO  " DAT  A 1  " 

200  BUFFER  #1 

210  FOR  1=1  TO  2 

220  LET  R<2, I )=R< 1,1) 

230  NEXT  I 

240  LET  J=2 

250  GOSUB  660 

260  LET  S=S+Q8 

270  LET  J=3 

280  FOR  1=1  TO  2 

290  LET  R ( 2 , I ) =  R ( 2 , I ) * G 6 

300  LET  R(3, I )=R( 1 , I )+R(2, I )/2 

310  NEXT  I 

320  GOSUB  660 

330  LET  J  =  4 

340  FOR  1=1  TO  2 

350  LET  R ( 3 , I ) =  R  < 3 , I ) *  G 6 

360  LET  R<4,  I  )=R<  1 ,  I  >+R<3,  I  >/'2 

370  NEXT  I 

380  GOSUB  660  •  1 

390  LET  J  =  5 

400  LET  S=S+Q8 

410  FOR  1=1  TO  2 

420  LET  R < 4 , I ) =  R < 4 , I)*G6 

430  LET  R(5, I )=R< 1 , I )+R(4, I ) 

440  NEXT  I 

450  GOSUB  660 

460  FOR  1=1  TO  2 

470  LET  R(5, I >=R(5, I )*G6 

480  LET  R  <  1 ,  I  )  =R  (  1  ,  I  )  h-  <  R  (  2 ,  I  )  +R  (  3 ,  I  >  *2+2*R<  4  ,  I )  +  R  (  5 ,  I  )  )  /6 
490  LET  R < 6 ,  I)=R<1,  I) 

500  NEXT  I 

510  PRINT  USING  100; S, 2, Si , T, P, Rho, Ust 

520  IF  Ust=0  THEN  210 

5  3  0  K  o  u  n  t,  =  K  o  u  n  t  + 1 

540  PR  I  NT  #  1  ; Ust , 1 , Rhost , Tst 

545  IF  Kount  MOD  400  THEN  210 

550  IF  NOT  Break  OR  (Rho >.8035)  THEN  210 

570  IF  SI  THEN  SI  a 

580  Sl=l 

598  PRINT  #1; 1E51, 0,8,0 
60O  G6=4"*G6 

605  08=4*08 

610  GOTO  210 

615  SI  a:  IF  Rho >.081  THEN  210 

628  PRINT  #1; 1E51, 0,0,0 

630  PRINT  # 2 J  0 , 0 , 0 , G 6 / 4 , O , A d a , 1, 1 . 25*SQR( Ast ) , 1 
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1329, .36? 


I 


648  L  0  H  Ei  "  n  e  w  n  s  n 11 
658  STOP 
668  F=R<Jfl) 

678  Si"RCJ,2> 

680  Z” ( l-F>*<C0+C2p*F> 


690  IF  Z>1  THEN  2=1 

780  IF  F >0  THEN  Z»< 1-F>*<C0tC2*F> 

718  IF  F < - . 5  THEN  2*1 
728  Z  h  e  1  p*  V0-C*  <  1  -2  >  s  Rh  op 
730  V  p  res*Zhe!  p  -  <  C  *  Z  3  i  >  *  R  d  a 
740  V  t.  h  *  Z  h  e  1  p -  C  *  Z  *  ft  d  a 

7  0  T  T  O  *  (.  V  t  h  *  <  2  -  S  i  /  C  >  /  <  V  p  r  e  s  *  Z  )  >  A  <  G  am  m  a~  1  > 

760  P  =  F  o  r  c  e  *  C  C  *  Z  *•  S  i  y  ^  T  f  <C  T  0  ■*  Vp  r  e «.  ) 

778  B r e a k  =  B r e  ale  OR  <  P  > P & t,  ar t  > 

788  R h o  =  <  C * Z - S  i  > /She  Ip 


7  9  0  S  i  gma~  R  h  o /  <C  1  -  R  d  aw  R  h  q  y 


308  R(.J,  1  >=-Bet  a*CP+Pt  > 

810  IF  NOT  Break  THEN  RETURN 

828  He  1 p  e  r  =  E  p  *  R  da*  3  i gma 

830  S  i  gma$t  -Set*  a*S  i  gma*  C  1  +  He )  per  > 

340  Rhost=Si qm$s%s< 1 +flda*S i gmast > 

8  5  w  T  s  x  -  2  *  T  *  <  1  +  C  G  am  m  a  - 1  >  *  H  *1  p  e  r  >  s  <  G  am  m  a+ 1  ) 
t* 0  P s t - P h el  p * P * (,  1  + 1 j am m a^ He  1  p e p ) 

870  U  si*  S  Q  R  <  3  8  6  *  G  a  m  fo  a  *  P  s  t  *  u  +  R  d  a  *  R  h  o  $  t  >  -  -  2  /  R  h  o  5.  t  ) 
880  R <  J ,  2 >  *fls t  *U=.t*Rhost 
890  RETURN 
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"newnsn" 


The  data  listed  in  this  program  should  not  normally  be  changed. 
Input  parameters  are  requested  as  indicated  below.  Default  values 
are  as  indicated. 


Message 


Keyboard  Entry  . 


0  FOR  GRAPHS,  1  FOR  NO  GRAPHS 
TIMM,  TEMPM 

ENTER  LABEL 
GUN  TEMPERATURE? 

ALCP 

XINT,  TINT 


NRDS,  NLA,  IPROPT,  LHELP 


INSIDE  DIAMETER, 
OUTSIDE  DIAMETER 


0  or  1 

Maximum  time  (sec)  and  temperature 
(°K)  to  be  shown  on  graph 
An  alphanumeric  label  for  the  graph 
Temperature  of  gun  steel 
(default  is  300° K) 

(default  is  0.4/(14.8  +  4*logio  D)^, 

D  is  the  gun  caliber  in  centimeters) 
Distance  from  idealized  breech  (inches) 
for  these  calculations,  time  to  stop 
calculations  or  go  to  next  round  (sec) 
(default  =  projectile  base  or  nozzle 
throat,  D/100) 

Number  of  rounds  in  this  burst,  number 
of  mathematical  layers,  print  option 
(0  for  full  print,  1  for  1/4  as  much, 

2  for  very  little  printing) ,  mathematical 
layer  to  use  for  bore  surface  between 
rounds  (default:  1,  7,  0,  4) 

Inside  and  outside  diameter  of  the 
gun  (inches) 
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1  OPTION  BASF  1 

5  DIM  Label *C 48] 

6  I  N  p  U  T  "  0  F  0  R  G  R  R  P  H  3 ,  ]  F  0  R  H  0  G  R  A  P  H  S  “  ,  E  r  i  q 

7  IF  Eng  THEN  355 

10  PLOTTER  IS  #l GRAPHICS" 

11  PEN  2 

12  DEG 

2  9  I  M  PUT  "  T I M  M ,  T  E  M  P  N  ",Ti  m  m ,  T  *  m  p  m 

2b  L INPUT  "ENTER  LABEL ", Labe  1  $ 

2b  Label  =LEN (Label  f  > 

30  LOCATE 

4  0  SCALE  -  T  i  m  m  1 0 ,  T  i  m  m  f  2  0  0 ,  T  e  m  p  m  +10  0 

45  CL  IR.  -T  i  mm/ 1 00 ,  T  i  mm  ,  290 ,  T empm 

5 0  A  X E  S  ,001, 1 0 0 , 0 , 3  © 0 

55  UNCLIP 

60  C SIZE  3 

70  LDIR  0 

30  LORG  5 

30  FOR  I = . 002  TQ  T i mm- , OO i  STEP  , 002 

100  MOVE  1,278 

110  LABEL  USING  " DD " ; I  * \ 088 

1 20  NEXT  I 

121  LDIR  98 
130  LORG  5 

1 4  8  F  0  R  1=4  0 0  T  0  T  e  m  p m  S T  E  P  200 
145  MOVE  83Q*Ti mm, I 
150  LABEL  USING  M B p D H “ ; J 
170  NEXT  I 

1  o  0  L  0  P  G  5 

1 9  8  M  0  V  E  -  ,  O  £  *  T  i  m  rn ,  (  3  8  8  +  T  e  m  p  m )  /  2 
2 0 0  C  SIZE  4 
281  LDIR  98 

210  LABEL  USING  " 28A" ; M TEMPERATURE  (DEG.  K > " 

228  LDIR  0 

2  3  8  M  0  V  E  T  i  m  m  /  2 , 2  3  8 

248  LABEL  USING  "  1  9 A  "  ;  ”  T  I  ME  (MILLISECONDS)  11 

2  5  0  M  0  V  E  T  i  m  m  /  2 ,  T  e  m  p  m  +  5  8 
268  C SIZE  6 , . 5 

278  LDIR  8 
298  LORG  5 
295  PEN  1 

308  LABEL  USING  “ 30 A " J Labe  1  $ 

310  PEN  3 
328  FRAME 
330  PEN  4 
340  MOVE  0,306 
355  OVERLAP 

361  D I M  Ac (9 9 ) , A  1 ( 9 9 ) 

37G  D  I M  I  st  he  <  7  >  ,  Th  <  99 ,4,7),  Hr  <  99 , 7  >  ,  I  st  a<  7),  I  o 

,Ak(40),Hold<7),Kno<7> 

371  DIM  Ah (4) 

372  MAT  READ  Ah 

373  DATA  . 25, 1 , . 5, . 333333333333 

3  7  4  I  p  r  o  p  t  =  1 

3  S  O  READ  T  h  t ,  E  o  p  r ,  L  h  e  1  p ,  T  m  a  x ,  \\  o  r  d  f 
38 1  MAT  Hoi d=<Tht > 

390  DATA  308, 1 , 4 , 6, 1 
40G  ASSIGN  #8  TO  "DATA" 

485  BUFFER  #8 

410  ASSIGN  #7  TO  " DATA 1 " 

411  Eof=0 

428  PE  AD  #7 ;  Rho  i  ,  wi  ,  T  i  ,  Du ,  T  ernp  i  ,  Ada,  Xm,  D ,  X  i  nt 
421  Xi=Xint 

438  INPUT  "GUN  TEMPERATURE? " , Tht 
440  MAT  Th  =  < Tht ) 

510  MAT  READ  I st a, I st o , ,?p , I st he 


J  p ( 7 ) , D  t( 7  > , D  a ( 7  > , H 1 ph<40) 
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DAT ft  1,5, 7,8,9, 9, 9,7, 11, 13, 15, 15, 15, 15,4,4,4,4,4,4,4,  1,5,7, 

y ,  9 , 9 

,9 

530 

REHD  Ere ad , Er , ErQ , Enew , Eos 

540 

DfiTfl  1 , 1 , 1 , 0,  1 

550 

Erh=0 

•560 

H2-1E50 

570 

E  r  h  =  0 

530 

R E H D  S,lrds,N r d s , K o s h 

599 

D  FIT  ft  0,  1  ,  1 , 0 

600 

Hi  c  p  ~  .  4  /  <  14.8+4*LG  T  ( D  >  >  A  2 

610 

INPUT  11  ft L CP? 11 ,  ft  1  c p 

620 

INPUT  " X I NT , T I  NT? " , X i  n  t ,  T  i n t 

6  3  0 

INPUT  “  NRDS ,  NLh  ,  I  PROPT  ,  LHELP  ? ,l ,  Nrds ,  N 1  a,  I  propt ,  Lhe  1  p 

640 

T  =  C!=Tc  al  =0 

/ 

s 

730 

L ini  =  L rn i  n  =  1 

740 

T int  »D*100 

' 

750 

N 1  a-7 

760 

MAT  READ  ftk  ,  ft  1 ph 

770 

DATA  6,0,0, 0,0,  •  108,  .  108,  .  187,  .  186,  .  105,  •  103,  .  101  ,  .098,  .096 

,  .095 

, . 083, . 033 

.  073 

, .073, .068, .062, .053, .06, .063, .089, . 187, . 122, . 130, . 135, . 140, 

.  145, 

.  145 

780 

D  FIT  ft  0, 0, 0, 0, 0, 0, 0, 0 

790 

DFlTft  0,0, 0,6,0,  .  128,  .  122,  .  115,  .  109,  .  103,  .096,  .09,  .084,  .078, 

.072, 

•  066 , • 06 1 , 

055, 

.049, .046, .045, .048, .056, .07, .089, . 1 13, . 128, . 135, . 140, , 145, . 

150,  . 

150 

800 

DftTft  0,0, 0,0, 0,0, 0,0 

820 

ft  1  phn>=0 

830 

FOR  K 1 = 1  TO  40 

840 

A1 phmsMAX<ftl phto, fll ph<Kl > ) 

850 

NEXT  K 1 

8  6  0 

Tdes i r  ®T i nt 

87G 

Rlce-422. 5*ftl  cp 

8  8  0 

Hspace=SQR<Du*ft 1 phm*3) 

890 

Const -Du/<Hspac e*Hspa»:  e ) 

900 

FOR  Kl«l  TO  40 

910 

R1ph(Kl > -Const*  F1 1 phCKl > 

920 

ftk  <K1  )  =  Const  *flk  <K1  > 

930 

NEXT  K1 

940 

Rlc  p*ft 1 c  e*Hspac  e*Const 

958 

PRINT  "THE  SPACE  I NTERVAL  IS",Hsp*ce 

960 

INPUT  "INPUT  INSIDE  DIftM,  OUTSIDE  DI AM" , Ri , Ro 

970 

Ri =2. 54*R  i 

980 

Ro=2. 54*Ro 

990 

K  o  u  =  8  -  N 1  a 

1 0  0  0  S  h  =  D  u  /  H  s  p  ac  e 
1010  D  t <  K  o u >  = D  u *  1 0 0 0 
1020  Ds (Kou ) =H space 
1038  FOR  Kol-Kou  TO  6 
1040  Dt  OCol  +  1 >=4*Dt <Kol > 

1 050  Ds < Ko 1  + 1 >  =2*Ds  < Ko 1 ) 

1060  NEXT  Ko 1 

1070  I s t o < 7  >  =  I  NT  <  < Ro-R i > /Ds  C 7  >  + .  5 > 
1080  REM  —CIRCULAR  CORRECTIONS 
1890  FOR  I  =  K o u  TO  7 
1100  1st Op= I S t  G  <1 ) 

111G  FOR  I i=l  TO  1st  op 
1120  R=R i+Ds<I>*Ii 

1130  DrCIi  ,  I>=Ds<I)/<2*R> 

1131  NEXT  Ii 

1132  NEXT  I 
1140  Eoe=  1 

1145  ft  1  ce=fl Ice*/ 12 
1150  FOR  K 7  = 1  TO  4 
1160  L=7 

1170  Kno < 7 >  =  J  =  K7 
1180  GOSUB  Heater- 
1190  IF  Eoe  THEN  J7 
1200  FOR  K  6  =  1  TO  4 
121G  L=6 
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1220  K  n  o  ( 6  )  =  J  =  K  6 
1230  G  O  S  U  B  H  0  a  t  0  r 
1240  IF  Eoe  THEN  J6 
1250  FOR  K5 - 1  TO  4 
1260  L  =  5 

1  2  7  0  K  n  o  5  )  —  J  —  K  5 
1280  GOSUB  Heater 
1290  IF  Eoe  THEN  J 5 
1300  FOR  K4=l  TO  4 
1310  L~4 

1320  Kno( 4  >  =  J  =  K4 
1330  GOSUB  Heater- 

1349  IF  Eoe  THEN  J4 

1350  FOR  K3--1  TO  4 
1360  L  =  3 

1370  Kno <3 ) =J  =  K3 
1380  GO S U B  H e at e r 
1398  IF  Eoe  THEN  J3 
1 4 O 0  F 0 R  K 2=1  TO  4 
1410  L  =  2 

1420  Kno < 2 >  = J  =  K2 
143  0  G  0 o U  £  H  e a t e r 
1440  IF  Eoe  THEN  J 2 
1450  FOR  K 1  =  1  TO  4 
1460  L= 1 

1478  Kno < 1 >  = J  =  K 1 
1 4  8  0  G  0  8  U  B  H  e  a  t  e  r 
1490  NEXT  K 1 
150Q  J2:  NEXT  K2 
1510  J3 : NEXT  K 3 
1528  J4 : NEXT  K4 
1530  J5 : NEXT  K5 
1540  J6: NEXT  K6 
155*0  J7 :  NEXT  K7 

1569  GOTO  1150 

1570  Heater: IF  Kosh< >0  THEN  2690 

1598  DEF  FNDe  1  (  I  a)  =Dr  <  I  a,  L  >  ~  <  fik  < I s 1 > -fik ( I s3 >  > / <  4*fik  (  I  *2  >  >  ; 

1  600  Ii  EF  FNF  (  fi ,  B ,  C  )=fi!(  K  >  *  <fi*  <  1  -Be  >-B-B+C*<  1  tBe  >  )+B 
1610  DEF  FNG (fi,  B,  C , E > =E*FNF < fi , B ,(!:>  +  <  1~E)*B 
1615  E  o p r = < I p r  o p t  =  1 )  0  R  ( I p r o  p t  =  2  >  fi  N  D  ( J  =  1 > 

1620  Jp<  L)= J 
1  6  3  0  R  u  =  0 

1640  Is  t  ar  t  =  I s t a ( L ) 

1650  Is t o p = I s t o ( L > 

1660  Il=Istart+l 

167G  12= I st op- 1 

1680  Jhe= 1  +  <  J  +  2)  MOD  4 

1690  IF  J=1  THEN  Hoi d(L>=Th< 1st op, 4, L) 

1700  FOR  1  =  I st art  TO  I  stop 
171G  L 1 =  T  h  < I , J  h e  ?  L  > / 5  0 
172G  F*  1  =  L  1  MOD  1 
1725  LI =1  NT (LI) 

1 730  fi 1 < I >  =fi 1 ph < L 1 > -P 1  * < fi 1 ph ( L 1 ) “R 1 ph < L 1 + 1 ) > 

1740  flc ( I >  =  fik (LI >-Pl *(fik (LI >-fik<Ll +1 > > 

1750  NEXT  I 

1766  IF  L  =  7  THEN  Out  side- 
1770  REM  — OUTER  STEP  OF  LfiY^R 
1780  I s= IstaCL+1) 

179Q  LI =Th( Is, Jp(L+l > , L+l )/50 
1  8  0  0  K  =  I  s  t  o  p 

1810  De-Dr  (K,  L)-(flc  (I  stop-2) -fik  <01  >-Pl  *(  fik  it  LI  >-fik  (LI  +1  )  >  >/<4*fic  <  1st  op)  ) 

1 820  Th ( I st  op , J , L >=FNG ( T h ( I st op-2 , 4 , L ) } Th< I st op , 4, L) , Th( Is, Jp<  L+ 1 ) , L+ 1 ) , fih ( J ) 
1830  Out  side: FOR  K=I1  TO  12 

1  y  4  u  D  e  -  D  r  (  K ,  L )  -  ( fl  c  (  K  - 1  )  —  fi  c  ( K  +  i  )  )  f  ( 4  *  fi  c  ( K  )  ) 

1  858  Th  (  K  ,  J ,  L  >  =FNF  (  Th  <K-1,  -Jhe ,  L  )  ,  Th  ( K ,  Jhe  ,L>,  Th  <  K+l  ,  Jhe-  ,L>> 

I860  NEXT  K  39 


1878  IF  L  =  7  THEN  Th < I st op , J , 7)  =  <2*Th< I  * t op- 1 , J , 7) -Th< I st op-2 , J , 7 > /2  + 1 . 8 ) / 1 . 566 
1888  REM  --THIS  ASSUMES  K= . 1 1 , H= . 08867 , AND  AMBIENT  TEMP  IS  300  K 
1898  IF  L  =  L Hi i  n  THEN  A48 
1988  I  s  =  I  st  o<L- 1 ) 

1918  Ll=ThOs-l,  l,L-l>/50 
1928  P 1 =L 1  MOD  1 

1938  De=Dr<  I  start ,  L>-<fik  tLl >-Pl  *<flk  <L1  >-Ak  <L1  +1  >  )-Ac  (  Istart+1  >  )/<4*Ac  <  1st  art  >  > 
1948  K= I start 

1 958  Th  < I st  art , J , L ) =FNF ( Ho  1 d< L- 1) , Th< 1st art , Jhe ,L),  Th  < I st  art  + 1 ,  The  ,  L  > ) 

I960  IF  Eopr  THEN  PRINT  USING  1 9©5;  S,  Rv, Q, L, ThC I st art , J, L) , Th< 1st  art  + 1 , J , L > , ThC I 
s  t  ar  t + 2 ,  J ,  L  >  ,  T  h  C  I  s  t  ar  t  +  3 ,  J ,  L  ) 

1965  IMAGE  DDD2 . D , 8D2 , 6DZ . DD , 4D , 7D2 . D, 7DZ. D, 7B2 . D , 7D2. B 

1978  GOTO  Helper  _  .... 

1988  A40: S=S+Dt <L) 

1998  IF  Eread  THEN  A58  .  " 

2880  IF  Eos  AND  CThC  1 ,  J,  LX  1 . 5*ThC5,  J,  L>  >  THEN  R12@ 

2810  R45:  ThC  1,  J,  L>=C2*ThC2,  J,  L>-ThC3,  J,  L>/2>/l . 5 
2020  Enew= C S >=Tdes i r ) 

2.030  IF  Eopr  THEN  PRINT  USING  1 965 ;  S ,  Ro ,  Sq ,  L  ,  Th  <  I  st  art ,  J  ,  L  >  ,  Th  <  I  st  ar  t  +  1 ,  J ,  L ) ,  Th  < 
Istart+2, J, L>, ThC Istart+3, J,L> 

2040  GOTO  Helper 
2845  A50:  M=M+1 

2046  READ  #8 ; V , X , Rho , Tgas 
2847  X=X+X i 

2050  IF  V  > 1 E  4  9  THEN  A78 
2860  A68: Vgas=0 

2078  IF  X>Xint  THEN  Vgas=V*Xi  nt  /X 
cd 8 8 8  H 1  =R  1  c  e*Rho*Vgas 
2898  Z=H1*H2/CH1+H2> 

2108  Coef  =  Z*Sh/Rc C 1  > 

2118  Th<  1 ,  J,  L)=C2*ThC2,  J,  L>-ThC3,  J,  L>/2+Coef *Tgas)/C  1 . 5  +  Coef  ) 

2111  Rv=Z*C Tgas-ThC  1 ,  J,  L  )  ) 

2112  IF  Eng  THEN  2128 

2113  IF  S>T i mm* 1000  THEN  Eng=l 

2119  DRAW  S/1000, Th O start , j, L> 

2128  Q =Q+CRv+R h o v > * D t CL>/2000 

2130  IF  Eopr  THEN  PRINT  USING  1 965; S, Rv , Q, L, Th< I st art , J , L ) , Th< 1st  art  + 1 , J ,  L >  ,  Th< I 


s t  a r  t  +  2 ,  J  ,  L  >  ,  T  h  Cist  ar- 1  + 

3,J,L> 

2140 

R  h  o  v  —  R  y 

2150 

GOTO  H  e 1 p e r 

2160 

R70: IF  E r  THEN  R80 

2170 

Er ead=0 

2180 

GOTO  ft 60 

2190 

fl80 : Er=0 

2200 

REM - STEP  SIZE 

INCREASE 

2210 

Sh-2*Sh 

2220 

fllO0:S=S-DtCL) 

2230 

Lrn  i  n=L  i  m=Mrn=L+ 1 

2240 

I i = I st  aCMm )  - 1 

2250 

•J.j  =  K  n  o  C  M  rrt ) 

2260 

FOR  1=1  TO  Ii 

2270 

ThC I , Jj,Mm>=ThC2*I 

-1,4, Mm- 1 > 

2275 

NEXT  I 

2280 

I  s  t  a.  C  M  rn )  =  1 

2290 

Kosh=3 

2300 

Kount  =0 

2310 

GOTO  2360 

2320 

H120:  IF  J<  >1  THEN 

H  4  5 

2330 

Kount  =  Kount  +  1 

2340 

EO'S=  <L<Lhe  1  p> 

2350 

ON  Kount  GOTO  R45, 

fl45, filOO 

2360 

He  1  per : Eoe=  <  Lm i n  >= 

L  > 

2370 

I  st  art  =  1st  aC  L  ) 

2380 

I  s  t  o  p=Isto<L) 

2390 

T  n*i  a x  =  M  R X  (  T  rri  ax  ,  T  h  (  1 

,  J ,  1  >  > 
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2408  IF  Enew=0  THEN  RETURN 

24  1  8  P’R  I  NT  "  Round  no"  ;  I  rds ;  “  T  I  HE=  " ;  3/  1 808  ;  "  Seconds " ;  UIH<  1  >  5  N  ax  t  f  nip "  ;  I  NT  < T max  > 
"Bore  surface  t  eniperat  ure  "  ;  Th<  i  ,  ,T?  Uhe  1  p>  ;  “Total  heat  i  nput  "  ;  ,  1*INT<  18*Q> 

2428  Q  =  0 
2430  Enew=0 

2448  IF  Nrds<=Irds  THEN  CALL  Uai  t <Hrd*. Tdesi r, Ti nt , Return) 

2450  IF  Return  THEN  RETURN 

2468  Tmax=Q 

2478  E o s  =  O 

2430  I  r d s -- 1  r  ds  +  1 

2490  Er ead=Er=Ero=True 

2500  Tdesi  r  =  Tdesi  r+T  i  nt 

2518  Sh“Sh/2 

2528  Lh=,L-l 

2530  FOR  1=1  TO  Lh 

2  5  4  8  J  p  <  I  >  =  4 

2558  NEXT  I 

2568  I=J 

2578  M m rn  =  L 

2575  E30 : Lh=M- 1 

2530  I  s  t  a  <  M  rn  rri  >  =  1  s  t  h  e  <  M  m  m  > 

2  598  I  s  =  I  s  t  a  »>  M  m  rn  > 

2608  FOR  Ka= 1  TO  Is  * 

2618  ThC2*K-l ,  4,  Lh)=Th<2*K,  1 ,  Uh>  =  <Th<K,  I ,  M  >+Th<K+l ,  I ,  M> ) / £ 

2628  Is t o p  =  I  s t o < L h ) 

2630  Ho  1 d <  Lh ) =Th < I s t op , 4 , Lh > 

2648  Mmm=Lh 
2650  1=4 

2660  IF  Lh>=2  THEN  E38 
2678  Lniin=Lim=l 
2630  RETURN 
2  6  9  0  K o s h = K o  s h- 1 
2788  E o e  = 1 
2718  RETURN 

2  7  2  0  S  U  E  14  a  i  t  •'!  N  r  d  s ,  T  d  e  S  i  r  >  T  i  n  t ,  R  e  t  u  r  n  > 

2  7  3  0  R  E  ft  Ii  N  a  dd ,  T  a  d  d 
2740  DATA  0,0 

2750  IF  Nadd=0  THEN  B10  ; 

2760  N r d  s = H a d d + N rds 
2770  Tint=Tadd 
2780  Return=0 
2790  SUBEXIT 

2300  B 1 0 : IF  Tadd=0  THEN  STOP 
23 10  T  des i r=Tdes i r+T  add 
2820  Return55! 

2830  SUBEND 
2840  END 
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3  Commander 
US  Army  Armament  Materiel 

Readiness  Command 
ATTN :  DRSAR-ASR 
DRSAR-LEA 
DRSAR-QAL 

Rock  Island,  IL  61299 
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DISTRIBUTION  LIST 


No .  of 

Copies  Organization 

1  Director 

US  Army  ARRADCOM 
Benet  Weapons  Lab 
ATTN :  DRDAR-LCB-TL 
Watervliet,  NY  12189 

5  Commander 

US  Army  ARRADCOM 
Benet  Weapons  Laboratory 
ATTN :  I .  Ahmad 

T.  Davidson 
G.  Friar 
P.  Greco 
J.  Zweig 

Watervliet,  NY  12189 

5  Commander 

US  Army  ARRADCOM 
Benet  Weapons  Laboratory 
ATTN:  J.  Busuttil 
W.  Austin 
R.  Montgomery 
R.  Billington 
J.  Santini 

Watervliet,  NY  12189 
1  Commander 

US  Army  Aviation  Research  & 
Development  Command 
ATTN :  DRSAV-E 
P.0.  Box  209 
St,  Louis,  MO  63166 

1  Director 

US  Army  Air  Mobility  Research  4 
and  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035 

1  Commander 

US  Army  Research  &  Technology  2 
Laboratories 
ATTN:  R.A.  Langsworthy 
Ft.  Eustis,  VA  23604 


No.  of 

Copies  Organization 

1  Commander 

US  Army  Communications  Rsch 
and  Development  Command 
ATTN:  DRDCO-PPA-SA 
Ft.  Monmouth,  NJ  07703 

1  Commander 

US  Army  Electronics  Research 
&  Development  Command 
Technical  Support  Activity 
ATTN:  DELSD-L 
Ft.  Monmouth,  NJ  07703 

1  Commander 

US  Army  Missile  Command 

ATTN:  DRSMI-R 

Redstone  Arsenal,  AL  35809 

1  Commander 

US  Army  Missile  Command 

ATTN:  DRSMI-YDL 

Redstone  Arsenal,  AL  35809 

1  Commander 

US  Army  Tank  Automotive  Rsch 
and  Development  Command 
ATTN:  DRDTA-UL 
Warren,  MI  48090 

1  Project  Manager,  M60  Tanks 

US  Army  Tank  &  Automotive  Cmd 
28150  Dequindre  Road 
Warren,  MI  48090 

Project  Manager 
Cannon  Artillery  Wpns  Systems 
ATTN:  DRCPM-CAWS 

US  Army  ARRADCOM 
Dover,  NJ  07801 

Project  Manager  -  M110E2 
ATTN:  J.  Turkeltaub 

S.  Smith 

Rock  Island,  IL  61299 
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No .  of 

Copies  Organization 

1  Project  Manager  -  XM1  Tank 
US  Army  Tank  Automotive 
Development  Command 
28150  Dequindre  Road 
Warren,  MI  48090 

1  Project  Manager 

Tank  Main  Armament 
ATTN:  A.  Albright 
Dover,  NJ  07801 

1  Project  Manager,  ARGADS 
Dover,  NJ  07801 

1  President 

US  Army  Armor  &  Engineer  Bd 
Ft.  Knox,  KY  40121 

1  Commander 

US  Army  DARCOM  Materiel 

Readiness  Support  Activity 
Lexington,  KY  40511 

2  Director 

US  Army  Materials  &  Mechanics 
Research  Center 
ATTN:  J.W.  Johnson 

K.  Shepard 

Watertown,  MA  02172 

3  Director 

US  Army  Research  Office 
ATTN:  P.  Parrish 

E.  Saibel 
D.  Squire 
P.0.  Box  12211 

Rsch  Triangle  Park,  NC  27709 

1  Director 

US  Army  TRAD0C  Systems 
Analysis  Activity 
ATTN:  ATAA-SL,  Tech  Lib 
White  Sands  Missile  Range, 

NM  88002 


No .  of 

Copies  Organization 

1  Commander 

US  Army  Air  Defense  Center 

ATTN:  ATSA-SM-L 

Ft.  Bliss,  TX  79916 

1  Commander 

US  Army  Armor  Center 
ATTN :  ATZK-XM1 

Ft.  Knox,  KY  40121 

1  Commander 

US  Army  Field  Artillery  School 
Ft.  Sill,  OK  73503 

5  Commander 

Naval  Surface  Wpns  Center 
ATTN:  M.  Shamblen 
J .  0 ! Brasky 
C .  Smith 
L.  Russell 
T.W.  Smith 

Dahlgren,  VA  22448 

2  Commander 

Naval  Ordnance  Station 
ATTN:  L.  Dickinson 

S .  Mitchell 

Indian  Head,  MD  20640 

1  Commander 

Naval  Ordnance  Station, 
Louisville 
ATTN:  F.  Blume 
Louisville,  KY  40202 

2  AFATL  (D.  Uhrig,  0.  Heiney) 

Eglin  AFB,  FL  32542 

1  National  Bureau  of  Standards 

Materials  Division 
ATTN:  A.W.  Ruff 
Washington,  DC  20234 
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No.  of 

Copies  Organization 

1  National  Science  Foundation 
Materials  Division 
Washington,  DC  20550 

1  Battelle  Columbus  Laboratory 
ATTN:  G.  Wolken 
Columbus,  OH  43201 

1  Lawrence  Livermore  Laboratory 
ATTN:  J.  Kury 

Livermore,  CA  94550 


2  Calspan  Corporation 
ATTN:  G.  Sterbutzel 

F.  Vassallo 
P.0.  Box  400 
Buffalo,  NY  14221 

1  Director 
Chemical  Propulsion  Info  Agency 
Johns  Hopkins  University 

ATTN:  T.  Christian 
Johns  Hopkins  Road 
Laurel,  MD  20810 

2  Princeton  University 

Forrestal  Campus  Library 
ATTN:  Tech  Lib 

B.  Royce 
P.0.  Box  710 
Princeton,  NJ  08540 

1  Purdue  University 

School  of  Mechanical  Engineering 
ATTN:  J.R.  Osborn 

W.  Lafayette,  IN  47909 


No.  of 

Copies  Organization 

1  SRI  International 

Materials  Research  Center 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

1  University  of  Illinois 
Dept  of  Aeronautics  and 
Aerospace  Engineering 
ATTN:  H.  Krier 
Urbana,  IL  61803 

Aberdeen  Proving  Ground 

Dir,  USAMTD 

ATTN:  H.  Graves,  Bldg.  400 

L.  Barnhardt,  Bldg.  400 
K.  Jones,  Bldg.  400 
R.  Moody,  Bldg.  525 

Cdr ,  TECOM 

ATTN:  DRSTE-FA 
DRSTE-AR 
DRSTE-AD 
DRSTE-TO-F 

Dir,  USAMSAA 
ATTN :  DRXSY-D 

DRXSY-MP ,  H.  Cohen 
D.  Barnhardt,  RAM  Div 
G.  Alexander,  RAM  Div 
Air  Warfare  Div 
Ground  Warfare  Div 
RAM  Division 

Dir,  USACSL,  Bldg.  E3516,  EA 
ATTN:  DRDAR-CLB-PA 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ _ _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name : 


Telephone  Number: 
Organization  Address: 


